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We present ab-initio electronic structure calculations for the Si-terminated SiC(0001)\/3 x \/H 
surface. While local density approximation (LDA) calculations predict a metallic ground state with 
a half-filled narrow band, Coulomb effects, included by the spin-polarized LDA+U method, result 
in a magnetic (Mott-Hubbard) insulator with a gap of f .5 eV, comparable with the experimental 
value of 2.0 eV. The calculated value of the inter-site exchange parameter, J — 30K, leads to the 
prediction of a paramagnetic Mott state, except at very low temperatures. The observed Si 2p 
surface core level doublet can naturally be explained as an on-site exchange splitting. 

PACS Numbers: 73.20.At, 75.30.Fv, 75.30.Pd, 7f.15.Mb 
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The fractional adlayer structures on semiconductors show 
rich phase diagrams with potential instabilities of charge 
density wave (CDW) and spin density wave (SDW) type 
at low temperatures, as well as Mott insulating phases 
§ . The (0001) V3 x ^3 surface of hexagonal SiC, (as well 
as the closely analogous (lll)v3 x \/3 of cubic SiC) ex- 
pected from standard LDA calculations to be a 2D metal 
with a half-filled narrow band of surface states in the 
bulk energy gap |^||] - is instead experimentally proven 
to be an insulator with a rather large (2.0 eV) band gap 
||,[|. It has been suggested that this system is a Mott 
insulator || due to the large value of ratio U/W, where 
U is the Coulomb interaction parameter, of the order of 
several eV, and W is the surface bandwidth, calculated 
to be about 0.35 eV P|| . Very recent STM data further 
confirm this picture |7||. However, a detailed electronic 
description of the ensuing state and of its magnetic im- 
plications is not yet available. 

Here we present a novel electronic structure calcula- 
tion for this system, based on the LSDA+U method ||, 
which is able to take into account Coulomb interactions 
between localized electrons. The result - an insulating 
surface, with a large intra-adatom exchange splitting, but 
an exceedingly weak antiferromagnetic inter-adatom ex- 
change coupling - is now in quantitative agreement with 
available photoemission data. 

The generalized LSDA+U method || is devised to par- 
tially cure the weakness of the local spin density approx- 
imation (LSDA) functional in dealing with interacting 
electrons in strongly localized orbitals, by supplement- 
ing the standard LSDA functional, .E LSDA [p <T (r)] where 
a =1,1, with a mean-field Hartree-Fock factorization, 
E u , of the screened Coulomb interaction V ee among elec- 
trons in localized atomic orbitals i, m. (Here i collectively 
labels the quantum number n, the angular momentum I, 
and the atomic site, while m labels the projection of the 
angular momentum.) The generalized LSDA+U func- 
tional is then written as: 
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\E u [{n^}]-E Ac , (1) 



= n^ m , is the localized orbital occupation 
density matrix. Double counting of the localized orbital 
contribution to the energy, already included by LSDA 
in an average way, is taken into account, in Eq. (]]]), by 
subtracting a compensating term 
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= ^W-1)-J^^-1)], (2) 



where U and J are the screened Coulomb and exchange 
parameters 0,0, N? = £ m ntf m , and N> = Nj + N % K 
Minimization of this modified functional with respect 
to the charge density p a (r) and the orbital occupation 



leads to an effective single particle Hamiltonian 



£ = £lsda + ]TE \ im <r)K. 
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which corrects the usual LSDA potential through a lo- 
calized orbital contribution of the form: 



C' = ^{(mV I Vee I m',m"')(n^, m „, 
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Important ingredients in the calculation are the matrix 
elements of screened Coulomb interaction T4 e which we 
parameterize, in analogy with the atomic case, in terms 
of effective Slater integrals F k |l2] which, in turn, can 
be linked to the Coulomb and Stoner parameters U and 
J, as obtained from LSDA-supercell constrained calcula- 
tions 11 13] . For Si, we apply the above corrections to 3p 



states only, since 3s states are fully occupied. The calcu- 
lated Coulomb and Stoner parameters for Si 3p states are 
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U = 8.6 eV and J = 1 eV, respectively. The electronic 
structure calculations we have performed rely on the Lin- 
ear Muffin-Tin Orbitals (LMTO) method @, and use 
the Stuttgart TBLMTO-47 code. An important point is 
that the calculation allows us, through the Green's func- 
tion method, to evaluate the effective inter-site exchange 
interaction parameters Jy, as second derivatives of the 
ground state energy with respect to the magnetic mo- 
ment rotation angle p5|,p|. 

The SiC(0001)-\/3x ^surface was simulated by a slab 
containing two SiC bilayers plus a Si adatom layer. The 
Si adatoms were placed in the T4 positions of the up- 
per (Si) surface. All atomic positions were fixed to the 
values calculated by Rcf. (2). The bottom surface (C 
atoms of fourth atomic layer) was saturated with hydro- 
gen atoms. A standard LDA calculation for the y/3 x V3 
unit cell containing one Si adatom gave, as expected, a 
metallic ground state with a half-filled narrow band sit- 
uated within the energy gap of bulk SiC (Fig. |l|). This 
surface band originates from the dangling bonds of the Si 
adatoms (p^-orbitals). The corresponding Wannier func- 
tion, however, contains less than 50% of the adatom p z 
orbitals, for, as is well-known it extends itself down 
into the first SiC bilayer. 

The LSDA+U method, relying on a Hartree-Fock-like 
mean-field approximation, can only deal with statically 
long-range ordered ground states. It is generally believed 
p6[ that a triangular lattice of electrons at commensurate 
filling, with nearest neighbor hopping, and a strong Hub- 
bard repulsion - i.e., in the Heisenberg limit - possesses, 
in spite of strong quantum fluctuations, a three-sublatticc 
120°-Neel long-range order - a commensurate spiral spin 
density wave with 120° spins lying on a plane. This kind 
of non-collinear magnetic order can in principle be han- 
dled by LSDA+U, through a straightforward extension 
of the functional. However, that extension makes the 
scheme not only computationally heavier, but of intrinsi- 
cally worse quality than the corresponding collinear cal- 
culation, as recently found in another case |L7| ■ We there- 
fore decided to carry out our calculation in the correct 
three-sublattice supercell (with three Si adatoms) but to 
restrict to collinear magnetic moments. Non-collinearity 
is then taken into account within a much less computa- 
tionally heavy Hubbard model, the resulting small value 
of antiferromagnetic coupling further supporting the su- 
perior accuracy of this approach. 

We started our LSDA+U calculation assuming two 
adatoms with a finite spin-up projection of the magnetic 
moment, and the third with spin-down. The result con- 
verges to a stable magnetic insulating ground state with 
an energy gap of 1.5 eV, which compares rather well with 
the experimental value of 2.0 eVJ9]- The corresponding 
energy bands are shown in Fig. pTOf the three spin-up 
bands, two are occupied and one is empty, (Fig. 0(a)), 
whereas for spin-down electrons (Fig. |^(b)) one band is 
occupied and two are empty. Consequently, there is a 



net magnetization corresponding to one spin for three 
adatoms. Of course, this net magnetic moment will dis- 
appear in the true noncollinear antiferromagnetic state, 
(except perhaps for a small fraction which might be sta- 
bilized by surface spin-orbit coupling). 

As a further check that the restriction to collinear spins 
does not really alter the main features of the bands in a 
substantial way, we performed a Hubbard model calcu- 
lation for the possible 3x3 magnetic phases, within a 
Hartree-Fock (HF) treatment [Q. Restricting ourselves 
to a one-band model with a nearest-neighbor hopping 
t w 0.04 eV (as extracted from the LDA surface band- 
width), we use an effective Hubbard U roughly equal to 
the ab-initio energy gap of « 1.5 eV. The results for the 
two magnetic solutions found are shown in Fig. |^(a,b). 
Fig. ||(a) shows the linear SDW solution with a uni- 
form magnetization m z — 1/3, which compares rather 
well with the ab-initio surface bands. Allowing for non- 
collinear spins, we then obtain the spiral SDW solution 
whose bands are shown in Fig. ||(b). We stress the fact 
that the spiral SDW is the actual HF ground state, while 
the linear SDW is only a metastable state. We note, how- 
ever, that the ground state energies of the two solutions 
differ by less than 0.5 meV/adatom, so that unrealis- 
tically precise ab-initio calculations would be necessary 
to decide which of the two is the actual ground state. 
Moreover, the HF bands for the two solutions are rather 
similar, apart from some extra splittings introduced by 
the collinear magnetic solution. It is not clear whether 
the present resolution of photoemission data would 
allow to distinguish between these two slightly different 
sets of bands. 

Coming back to the ab-initio study, we have calculated 
a magnetic moment for the Si adatom of approximately 
0.6/iB- This reduction is of course not due to quantum 
fluctuations, absent in the LSDA+U, but to the subsur- 
face derealization of the Wannier function, whose weight 
is only about half on the Si adatom. For the same rea- 
son, the energy gap value of 1.5 eV is roughly 4—5 times 
smaller than the bare U value. In fact, the LSDA+U 
correction to the potential of a particular orbital ma is 
roughly proportional to 17(1/2— n m(7 ), and the splitting of 
the potentials for spin-up and spin-down p z -orbitals of Si 
adatoms is proportional to the magnetic moment on the 
adatom. This means that this potential splitting can be 
estimated to be roughly U/2. If one takes into account 
that the potential correction is applied to Si3p-orbital, 
which is about one half of the total Wannier function 
corresponding to the half-filled band, its effect on the en- 
ergy splitting AE will be further reduced by a factor two, 
yielding roughly AE ss [7/4, close to the value obtained 
in our actual LSDA+U calculation. 

We have also calculated the value of the inter-site ex- 
change parameter corresponding to the interaction 
between the spins of the electrons localized on neigh- 
boring Si adatoms, and obtained a value of Jy = 30K. 
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The corresponding one-band Hubbard model estimate is 
Jij = 4t 2 /U e s, where t is the nearest- neighbor hopping 
parameter and U e s is the effective Coulomb Hubbard U. 
Taking U c s of the order of the energy gap (« 2 eV) and 
t = 0.04 eV, we obtain 37K, which agrees quite well with 
the LSDA+U result. 

The existence of finite magnetic moments on the Si 
adatoms should be experimentally detectable. Meth- 
ods which require magnetic long-range order do not 
seem viable, since the only circumstance where a finite- 
temperature order parameter could survive in this frus- 
trated 2D magnetic system, would be in presence of 
a large spin-orbit induced magnetic anisotropy favor- 
ing out-of-surface ordering, which we do not have rea- 
son to expect. Therefore, without ruling out the possi- 
bility of low-temperature ordering, we believe that the 
SiC(0001)V3 x V3 surface will be in an overall param- 
agnetic Mott insulating state, in spite of the existence 
of an on-site moment, at least down to liquid nitrogen 
temperature. 

The situation is much more promising at the intra- 
atomic level. Here, exchange splittings in the Si adatom 
2p core levels should be large and detectable. One com- 
plication is that the adatom dangling bond has a pre- 
dominant 3p z character, which breaks the core level 
symmetry between 2p z and 2p XiV by an amount not 
negligible in comparison to spin-orbit coupling and ex- 
change. The 2p core hole will exhibit a 6-fold multiplct 
resulting from the joint effect of intra-atomic exchange, 
asymmetry, and spin-orbit. An estimate for the multi- 
plet splitting is obtained by assuming the valence elec- 
trons to be distributed with weights a in 3p z |j an d 
(3/A in each 3p x . y fj.. For Aso = 0, we calculate a 
bare asymmetry splitting A = 62 Px -\ ~~ e 2?^T — a [(Jzz ~ 
Jzx) - {K zz - K zx )\ + {/3/4)[2(K zz - K zx ) - {J zz - J zx )}, 
where K zz(x) = (2p z{x) , 3p z |e 2 /|i"i - T 2 \\2p z{x ), 3p z ), 
and J zz(x) = (2p z(;E ),3p 2 |e 2 /|ri - r 2 \\3p z ,2p z(x) ) , are 
Coulomb and exchange integrals. The first square 
bracket dominates, and is estimated by a Si + Hartree- 
Fock calculation to be about 1 eV, so that A ss a (eV). 
Since only about 50% of the dangling bond orbital is 
3p z , we obtain the desired crude estimate A « 0.5 eV. 
Exchange splittings are large for the 2p z case and small 
for the 2px,yi reflecting the large differences in the cor- 
responding exchange integrals J zz and J zx . When we 
include the spin-orbit interaction AsoL ■ S, the final core 
hole levels are obtained as in Fig. § For Aso = 0.4 eV we 
find a roughly three-peaked structure, the broad central 
peak four times as strong as each side peak. This offers 
an alternative explanation of the experimental lineshape 
p8[ (see Fig. |J) in terms of a single exchange-split mul- 
tiplet, rather than two chemically inequivalent sites Si 
and S 2 |@ whose existence is otherwise not supported. 

Additional direct experimental evidence for the Mott- 
Hubbard state could be obtained by a careful study 



of adatom (3p z , 3p z ) Auger spectral intensities, which 
should be easily singled out owing to the large gaps. The 
probability of double occupancy of the adatom dangling 
bond orbital should be almost completely suppressed, 
dropping from the band value of 1/4, to a value of or- 
der t/U which is one order of magnitude smaller. Even 
considering that only half of the orbital is adatom 3p z , 
this surface should show a (3p z , 3p z ) Auger intensity 
which, by comparison with the remaining (3p,3p) val- 
ues, is anomalously small, as a proof of its Mott-Hubbard 
state. 

More experimental and theoretical effort is clearly 
called for to check these strong correlations and related 
magnetic effects, possibly the first of this magnitude to 
be suggested for an sp bonded, valence semiconductor 
surface. 
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FIG. 3. Tight-binding HF electronic bands along high 
symmetry directions of the Brillouin Zone for the one-band 
Hubbard model at U c s/t\ = 40 for the two magnetic so- 
lutions: (a) linear SDW with a net uniform magnetization 
m z = 1/3; (b) spiral SDW, with finite on-site magnetization, 
and m z = 0. Solid and dashed lines denote up and down 
bands, respectively. 



FIG. 1. Energy bands obtained from a standard LDA cal- 
culation for the SiC(0001)\/3 x v3 surface. Brillouin zone no- 
tations correspond to a unit cell with one Si adatom. Energy 
is measured from the Fermi level. Note the narrow half-filled 
surface band. 
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FIG. 2. LSDA+U energy bands for the spin-polarized 
collinear state of SiC(0001)3 x 3 surface. Brillouin zone nota- 
tions correspond to the unit cell with three Si adatoms. (a): 
spin-up electrons, (b): spin-down electrons. Energy is mea- 
sured from the Fermi level. Note the insulating state, with a 
gap of about 1.5 eV. 



FIG. 4. Si 2p core levels calculated as a function of 
spin-orbit coupling, in comparison with photoemission data 
and original fitting in terms of bulk (Bl) and two surface 
sites SI and S2 by Johansson et al. Jl8| . Upon inclusion of in- 
tra-atomic exchange splitting the whole surface contribution 
S1+S2 can be explained as due to a single site. 
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